Reducing multiphoton ionization in a linearly polarized microwave field by local 

control 
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We present a control procedure to reduce the stochastic ionization of hydrogen atom in a strong 
microwave field by adding to the original Hamiltonian a comparatively small control term which 
might consist of an additional set of microwave fields. This modification restores select invariant tori 
in the dynamics and prevents ionization. We demonstrate the procedure on the one-dimensional 
model of microwave ionization. 

PACS numbers: 32.80.Rm, 05.45.Gg 



I. INTRODUCTION 



The multiphoton ionization of hydrogen Rydberg 
atoms in a strong microwave field is an experi- 
ment which revolutionized the way we view the physics 
of highly excited atoms Q (for thorough reviews, see 
^ IE IE Q)' Its interpretation remained a puzzle un- 
til its stochastic, diffusive nature was uncovered through 
the then-new theory of chaos |E IE H H 113 ■ 

The multi- 
photon ionization is believed to occur when the electrons 
diffuse to increasingly higher energies chaotically by tak- 
ing advantage of the breakups of local invariant tori in 
phase space IE IE ^1 13,, 14,2^]. Classical theory can 
be applied to the ionization of hydrogen in the parame- 
ter regime where the microwave and Kepler frequencies 
are nearly equal 0, 0| • Because the classical dynamics 
of this system is chaotic, the Rydberg states of hydro- 
gen are an excellent testbed for investigatiiig the quan- 
tal manifestations of classical chao s lla. llTl llq . i.e., the 
field of "quantum chaology" |E IE Il5l Il9j|. Indeed, the 
literature on the correspondence between classical and 
quantum behavior in the ionization of Ry dberg atoms is 

extensive JlJE BiJE IIEIIilillIllIllIEl!ElillIE 

[19, 20, 211121 l2al2^. We perform our purely classical 
calculations in the regime where the quantum-classical 
correspondence is particularly close ItI . 

Recently, the research focus in this field has shifted 
from understa nding to manipulating the ionization pro- 
cess |2E|2^,'27, 28, 29, "so]. Since microwave ionization of 
Rydberg states is a paradigm for time-dependent nonin- 
tegrable systems, learning to manipulate stochastic ion- 
ization is expected to pave the way to controlling other, 
more involved systems. The control of stochastic ioniza- 
tion has been investigated using both quantum and clas- 
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sical approaches in the past few years |2fil . 1271 1281 l29ll3H . 
Here, we return to the basic dynamics of the stochastic 
ionization process to answer the most elementary manip- 
ulation question: If the multiphoton ionization is made 
possible by broken invariant tori, can ionization be re- 
duced (or even stopped) by restoring invariant tori at 
carefully chosen locations in phase space? 

In this paper, we will show how to reduce or shut off the 
ionization of Rydberg atoms using a "local" control strat- 
egy which originates in plasma physics |32j | . The premise 
of the procedure is to reduce the chaos (and thus ioniza- 
tion) in a selected parameter range through a small per- 
turbation which regularizes the dynamics in that narrow 
area but does not affect the dynamics elsewhere. Techni- 
cally, local control achieves this by creating an invariant 
torus in a selected region of phase space without signifi- 
cantly changing other parts of phase space. 

The problem of finding such a modification of the 
original Hamiltonian system is, a priori, nontrivial : A 
generic modification term would lead to the enhancement 
of the chaotic behavior (following the intuition given by 
Chirikov's criterion |^). Modification terms with a reg- 
ularizing effect are, of course, rare. However, there is a 
general strategy and an explicit algorithm to design such 
modifications which indeed drastically reduce chaos and 
its attendant diffusion by building barriers in phase space 
js^ls^ l. as we will show on the one-dimensional hydrogen 
Rydberg atom in a microwave field. 

One-dimensional models of microwave ionization in 
linearly polarized microwave fields have proven per- 
fectly adequate to explain most experimental observa- 
tions ^4, j IE since many of the experiments con- 
sidered extended, quasi-one-dimcnsional hydrogen atoms 
0, 0, 0i in which the angular momentum of the Ryd- 
berg electrons is much smaller than their principal quan- 
tum number. As a result, the atoms resemble needles in 
which the electron bombards the core with zero angular 
momentum. The Hamiltonian in atomic units reads 



H{p,x,t) 
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Xx cos tut, 
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where A is the ampUtude of the external field. 

The desired Hamiltonian with the control field reads 

H(p,x,t)^^ + Xxcosujt + xf(t). (2) 

2 X 

For practical purposes, we expect that the control field 
f{t) has the same form as the perturbation field but with 
relatively small amplitude. Despite the fact that / intro- 
duces an additional set of resonances, its effect, if it is 
appropriately chosen, is to restore specific invariant tori. 

This paper is organized as follows: In Sec. m after 
summarizing the control method [33 |. we implement it 
on a one-dimensional hydrogen atom driven by a lin- 
early polarized microwave field. In Sec. IIIII we present 
the numerics of the control term and show its efficiency 
by using Poincare sections, laminar plots and diffusion 
curves. In order to be relevant and feasible for physi- 
cal implementations, the control term has to be robust, 
i.e., sufficiently good approximations to it should reduce 
chaos effectively, too. We pay particular attention to this 
point and show numerically that reasonable approxima- 
tions to our control terms are effective in reducing chaos 
also. Conclusions are in Sec. IIVI 



II. COMPUTATION OF THE CONTROL TERM 

We consider Hamiltonian systems with L degrees of 
freedom, written into action-angle variables (A, 9) of the 
form 

H{A,9)^Ho{A)+eV{A,9), (3) 

where {A, 9) e x T^. 

In the integrable case (e = 0), the phase space is foli- 
ated by invariant tori with frequency '^'(A) = dHo/dA. 
Let us consider one of these invariant tori with fre- 
quency u) and position Aq. We assume that u) is non- 
resonant, i.e. there is no non-zero integer vector k such 
that w-k = 0. This invariant torus is generally destroyed 
by the perturbation T^(A, 9) when the parameter e is 
greater than a critical value Sc- The idea of the control 
is to rebuild this invariant torus with frequency lj for 
e > Ec by adding a small control term / to the original 
Hamiltonian H. Thus the controlled Hamiltonian He can 
be constructed as 

H,{A,9) = H{A,9) + f{e). 

The expression of control term is given by 

f{9) = -H{Ao~deTb{9),9), (4) 

where b{9) = H(Ao, 9) and F is a linear operator defined 
as a pseudo-inverse oi us ■ de. Its explicit expression is 

Tb{9)= V -^e^i^^ (5) 



for h{0) = X^fcGZ^ b^e^^'^. The restored invariant torus 
of the controlled Hamiltonian He has the equation : 

A = Ao-rdeH{Ao,e). (6) 

Such an invariant torus acts as a barrier to diffusion for 
Hamiltonian systems with two degrees of freedom. 

In order to apply this method to a one-dimensional 
hydrogen atom driven by a microwave field, we first need 
to map Hamiltonian |^ into action-angle variables of the 
unperturbed system (A — 0). Its action-angle variables 
(J, 9) are [13 

X = 2 sin^ 95, 
1 

p= -jCOt If, 

with 

9 = 2(p — sin 2ip. 

After rescaling energy, time, position and momentum as 
H' = uj-'^/^H,t' = ujt,x' = J^l'^x.v' = bj-'^/'^p, we 
obtain the rescaled field amplitude A' = uj~^/^X. The 
rescaled Hamiltonian still satisfies the equations of mo- 
tion, and we assume uj — 1 without loss of generality in 
Eq. 1^. The scaled frequency, or say, winding ratio in 
the rescaled system is thus defined as = J^u — . 
Expanding x', we rewrite Hamiltonian Q |^ 

^ 00 
i? = 2 j2A(y + ^a„ cos n6') cost, (7) 

n— 1 

where 

_ Jn{n) - Jn-ijn) 

— 1 

n 

and J„'s are Bessel functions of the first kind. We ab- 
breviate the Hamiltonian lO as 

H = -^ + \J^v(9,t). 

where 

00 

v{9, t) — ao cos t + a„ [cos (n6> + t) + cos{n9 — t)]. (8) 

n=l 

The Hamiltonian {Tj) displays a set of primary resonances 
approximately located at J„ = n^^^. The overlap of these 
resonances 33] leads to large-scale chaos and hence ion- 
ization. We expect the lower action region to be more 
regular, and it is chaotic for sufficiently large A by res- 
onance overlap. The idea is, given a value of n, to re- 
store an invariant torus in between the resonances ap- 
proximately located at J„ and Jn+i- For A = 0, this 
invariant torus with (Kepler) frequency luq is located at 

— 1/3 

Jq = lOq .In order to do this, we compute the control 
terms as explained above. 
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The next step of the control algorithm is to map the 
time-dependent Hamiltonian into an autonomous one. 
We consider that t (modulus 27r) is an additional an- 
gle variable and wc call its corresponding action variable 
E. The autonomous Hamiltonian becomes H{J, 9, t)+E. 
The action-angle variables are A — (J, E) and 9 = {9, t). 

The frequency vector of the torus is u; = (o^q, !)■ The 
formula of the control term is obtained by replacing the 
actions A by Ao-9eri?(Ao, 9) where Aq = (Jo, 0). The 
control term / is given by 



III. NUMERICAL ANALYSIS 

In what follows, the series which give v and Tdgv are 
truncated at n = 30 for numerical purposes, and the first 
series of / is truncated aX k = 20. Also, we choose ujq 
in the interval { ^^^^i , which corresponds to a region 
in between two primary resonances, and ng is in general 
chosen equal to 1, 2, 3... (With a relatively big no, local 
control theory still holds though quantum suppression 
leads to a higher ionization threshold "sB*!). 
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fc=2 



+ {2\^uj-^Tdgv - X^uj^ ^ {Vdevf)v, (9) 



where 



Tdev = ^ nan 



cos{nd + t) cos{nd — t) 
nujo + 1 nujo — 1 



For A small, we approximate the control term / by its 
leading order in which is given by 



f2{9,t) = -X\TdgVy + 2X'LU^'vTdgV. 



(10) 



Obviously, the location of the restored invariant torus 
depends on the choice of Kepler frequency loq or equiva- 
lently of its location in the integrable case Jo. The theo- 
retical torus curve is given by 



J = J{9,t) = ^0 ' - ^^0 'Tdev{9,t). 



— 1/3 

This torus is A-close to Jo = ij^o • We notice that the 
control term as well as the invariant torus are 27r-periodic 
in 9 and time t. 

Remark : In the local control method, we have searched 
for control terms only dependent on 9. However, in order 
to be more consistent with the specific shape of the con- 
trol waves, it can be appropriate to search for controlled 
Hamiltonian of the form 



(11) 



H,{A,9) = H{A,9) + {n-A)^f{9), 



where is a fixed vector, e.g., fi = (1, 0) 
Following the same arguments as in Ref. [32 
of the control term is 



fiO) 



g(Ao - dgTb, 9) 
{Vl-Ao-n- deThf 



in that case. 
I , the formula 



(12) 



where h{9) = H{Aq,9). We notice that the control 
term H12|l is still of the same order as the one given by 
Eq. jnj and it is £^-close to the one given by Eq. 
divided by (Jl • Aq)^, since dgb is of order e. 



A. Analysis of the control term 

Figure n depicts a contour plot of / given by Eq. 
and /2 given by Eq. (|10|l for luq = 0.6750 (which corre- 
sponds to no = 1) and A — 0.03. In this case, the scaled 
frequency at the intended invariant torus for A = is 
ru = Jq = ujq^ = 1.4815 < 2 which justifies the appli- 
cation of classical theory to the regime we are interested 
[Tl . 36] . Since / and /2 are 27r-periodic in t and 9, these 
contour plots are represented for {t,9) £ [0,27r]^. In or- 
der to compare the control term with the perturbation, 
Fig. El represents a contour plot of the perturbation at 
an action J = Jq where the control acts. These figures 
show that for this value of A the control term is small (by 
a factor approximately equal to 10) compared with the 
value of the external field XJQv{9,t). 

However, the control terms / and /2 given by Eqs. © 
and (|10|l appear to have a much richer Fourier spectrum. 
We have represented in Fig. |31 their two-dimensional 
Fourier transforms. They have an infinite number of 
Fourier modes and therefore not practical for a numerical 
or experimental realization. However, it is seen on Fig.|31 
that only few Fourier coefficients contribute significantly 
to the control terms. Therefore, it is feasible to truncate 
them since the method has been shown to be robust [s^l ■ 
The tailored control term results in general from a trade- 
off between the ability to control chaos and restrictions 
on the desired shape for the specific problem at hand. 

In order to identify the main Fourier modes, we intro- 
duce a parameter A defined as 

A _ l/fcl,fc2l 

^fel,fc2 = 



\kiUJo + k2\ 

where fki,k2 is the Fourier coefficient with wavevector 
(fci,A;2) of / or /2. The dominant Fourier mode is sup- 
posed to have maximal A. We notice that this definition 
contains two effects : First a dominant Fourier mode has 
to have a significant amplitude, and second, its corre- 
sponding wavevector has to be close to a resonance with 
the frequency vector of the integrable motion (and hence 
close to a resonance). For loq = 0.6750 (which corre- 
sponds to no = 1) and A = 0.03, there is only one domi- 
nant Fourier mode in / or /2 which has a frequency which 
is twice the microwave frequency. The truncated control 
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(a) 




(b) 



FIG. 1: Contour plots of (a) / given by Eq. © and (b) /2 given by Eq. (EU for A = 0.03 and cjq = 0.6750. 




FIG. 2: Contour plot of XJ^v{e,t) where v is given by Eq. Q for A = 0.03 and wo = 0.6750. 



■ 
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2 (a) 2 (b) 

FIG. 3: Two-dimensional Fourier transforms of (a) / given by Eq. @ and (b) /2 given by Eq. 1101 for A = 0.03 and luq = 0.6750. 



term is given by 

/a(^,t) =/3,-2COs(30-2t), 



(13) 



where /3,_2 ~ —9.772 x 10^* for control term / given by 
Eq. lO and /a, -2 ~ —9.739 x 10"'' for the approximate 
control term /2 given by Eq. Hl()|) . We notice that these 
two values are very close. For this mode, we have ^3,-2 ~ 
3.90 X 10~^ which is more than ten times larger than the 



second largest one A1.-2 ~ 1-90 x 10~^. We notice that 
the continued fraction expansion of loq is [0, 1,2,...]. One 
best approximant is [0, 1, 2] = 2/3 which is the frequency 
of the mode of fa- 
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FIG. 4: Poincare sections of (a) The uncontrolled Hamiltonian H given by Eq. 0, (6) The controlled Hamiltonian H + f 
where / is given by Eq. @, (c) The controlled Hamiltonian H + f2 where /2 is given by Eq. IIUII . and (d) The controlled 
Hamiltonian H + fa where fa is given by Eq. 11311 for A = 0.03 and uiq = 0.6750. The thin wary curve indicates the location 
where the invariant torus is restored. The black dots are from trajectories launched below this curve, and gray dots are from 
trajectories launched above this curve. Note how they are interspersed in (a), as is expected of chaotic trajectories, and how 
the control restricts their movements in phase space through the invariant torus. 



B. Poincare sections 

In order to test the efficiency of the control terms to 
restore invariant tori in phase space, we perform Poincare 
sections oi H + f, H + f2 and H + fa and compare them 
to the Poincare section of H given by Eq. ([TJ). Since all 
these Hamiltonians are periodic in time with period 27r, 
the natural Poincare section is a stroboscopic plot with 
period 27r. 

Figure 01 depicts Poincare sections of Hamiltonian {Tj) 
in panel (a), Hamiltonian H + f where / is given by 
Eq. lO in panel (6), Hamiltonian H + f2 where /2 is 
given by Eq. H10(l in panel (c) and Hamiltonian H + fa 
where fa is given by Eq. ifT^ in panel (d) for luq = 0.6750 
and A = 0.03. We notice that with the addition of the 
control terms, an invariant torus has been restored which 
prevent the diffusion from below to above the invariant 
torus. It is also worth noticing that all of these con- 
trol terms are efficient although only / is expected to 
be, indicating that the presence of the control field /2 
contributes dominantly to a restoration of invariant tori 
at specific locations such that higher order resonances 
are eliminated which, in the Chirikov's approach |33| . 
leads to less chaos and hence less stochastic ionization 



|32l l37j in our problem. It reinforces the robustness of 
the method and allows one to tailor a control term which 
is simpler to implement. 

C. A control term as an additional wave 

In Sec, nil A1 we show that the frequency of the control 
wave should be twice the one of the initial wave 
Therefore, a possible controlled Hamiltonian is 

p2 1 

H= h Ax cos i + /ix cos 2i, (14) 

2 X 

which corresponds to a control terms g{t) = ^cos2t in 
Eq. 121). In order to obtain the value of /i, we use the 
Fourier decomposition of the control term / obtained 
previously. First we map the controlled Hamiltonian into 
action-angle variables : 

H = + 2J2 {Xcost + ^icos2t)v{e,t), (15) 

where v is given by Eq. ijHJl . If we give a value of the action 
Jo where the invariant torus has to be restored, we have 
seen that the dominant Fourier mode is proportional to 
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FIG. 5: Poincare sections of Hamiltonian 11511 for A = 0.03 and ji = 0.0127. 




(a) 




(b) 



FIG. 6: Laminar plots of (a) Hamiltonian Q and (b) Hamiltonian l)15|l for A = 0.03 and /i = 0.0127. Cut-off time is GOOtt and 
diffusion threshold is Jth = 1.30. 




FIG. 7: Mean diffusion time (Td) versus initial action J with amplitude of external field A — 0.03. Solid lines are for 
Hamiltonian Dash-dotted lines are for Hamiltonian I15I I for jj, — 0.0127. Cut-off time is GOOvr and diffusion threshold is 
Jth = 1.30. 



cos[(2n+l)0 — 2i] where n is obtained using the continued 
fraction expansion of = •^o"'^- This mode is present in 
Eq. 1151) and has an amphtude given by ^jQa2n+i- If 
/2n+i,-2 denotes the amplitude of the dominant Fourier 
mode in Eq. for the values of the parameters luq and 
A, then the amplitude of the control field is chosen to be 



The real parameters taken for the external microwave 
field and control field are flexible for a set of rescaled 
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amplitude of external field A = 0.03 and frequency uj = 1 
as long as they satisfy rescaling relationships. 

Figure depicts the Poincare section of the controlled 
Hamiltonian ^ with A = 0.03 and ^ = 0.0127. FigureEl 
docs not show the restoration of an invariant torus as in 
the previous cases and hence, the ionization reduction is 
not obvious. This comes from the fact that the additional 
wave is quite far from the control term ij^ due to addi- 
tional resonances which break-up the restored invariant 
torus. 

However, the ionization process is still reduced, and 
this can be seen by looking at laminar plots. Such plots 
are obtained by looking at a grid of initial conditions and 
plotting the number of iterations it takes the action to 
exceed a certain threshold. Figure El depicts the laminar 
plots for Hamiltonian and Hamiltonian (|15|1 with A ~ 
0.03 and /i — 0.0127. The action threshold is chosen to 
be Jth — 1.30. The maximum integration time is GOOtt. 
The darker the region is the smaller time it takes to have 
J > Jth- It is expected that the laminar plots with more 
brighter regions are cases where there is less ionization. 

In order to compare the diffusion time of trajectories 
for Hamiltonian {T)) with that of the controlled Hamilto- 
nian H15|l , we have taken a set of JV initial angles evenly 
distributed in [0, 27r] for one initial action J and then 
computed the mean diffusion time for each J in both 
controlled and uncontrolled cases : 

1 ^ 

{Td){J)^J^J2^,{J,0^). (17) 

1=1 

Figure depicts the the curve of mean diffusion time 
{Td) versus initial action J. In the numerical computa- 
tion of Td{J, 6i), the integration is performed till the cut- 
off time t — GOOtt. Therefore for some trajectories the 
actual diffusion time is certainly above the cut-off time 



or even goes to infinity. The double frequency control 
field also works for the regime 1.14 < J < 1.20 for the 
reason that the rebuilt invariant torus is a curve which 
goes beyond Jq — ujq =1.14 for some areas. Figure [7| 
shows that the mean diffusion time for controlled Hamil- 
tonian H15|) is significantly larger than that for Hamilto- 
nian {Tj) which clearly shows the effect of the additional 
microwave field to reduce ionization. 



IV. CONCLUSION 

In this paper we implemented a local control method 
on the one-dimensional hydrogen atom in a linearly po- 
larized (LP) microwave field in order to reduce ionization. 
After simplifying the originally complicated control func- 
tion numerically, we obtained an extremely simple con- 
trol term which is in the same form as the external LP 
microwave field but with smaller amplitude. Adding the 
small control field to the perturbed Hamiltonian leads to 
a reduction of ionization. We have done the calculations 
in a regime where the quantum and classical agree, and 
our classical computations show efficient suppression of 
ionization. Preliminary results we obtained for control- 
ling ionization in higher dimensions show the promise of 
the local control method for manipulating ionization in 
full-dimensional Rydberg atoms and work in this direc- 
tion is currently under way in our center. 
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